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ABSTRACT 



The forces acting on yawed smooth and rough circular 
cylinders in oscillating flow have been investigated for the 
purpose of determining the appropriate force-transfer coeffi- 
cients and the applicability of the "independence principle.” 
The results have shown that the flow about each cylinder is 
unique and the independence principle does not hold true over 
the range of Reynolds numbers and Keulegan-Carpenter numbers 
covered by the investigation. It has further been shown that 
Morison's equation predicts the measured force with the same 
degree of accuracy as that for the normal cylinder provided 
that the force-transfer coefficients appropriate to each yaw 
angle, Reynolds number, Keulegan-Carpenter number, and the 
relative roughness are used. 



3 



PREFACE 



The work reported herein is part of a basic research 
on "yaw and current effects on the hydrodynamic resistance 
of smooth and rough cylinders in harmonic flow" and was 
performed during the period 1 September 1982 to 1 March 
1983. The principal investigator is Dr. T. Sarpkaya, Dis- 
tinguished Professor of Mechanical Engineering. The project 
was supported by the Water Resources and Environmental Engi- 
neering Program of the National Science Foundation under 
Agreement No. CEE-8210246. The cognizant NSF Program Officer 
is Dr. Arthur A. Ezra. 

The writer wishes to acknowledge the generous support 
of the National Science Foundation and the constructive comments 
of Dr. Arthur A. Ezra. A special note of thanks is extended 
to Mr. Thomas S. Raines and Mr. Dean 0. Trytten for their 
assistance with the experiments during the course of the 
investigation . 



4 



TABLE OF CONTENTS 



I. INIRCDUCTICN 12 

II. EXPERIMENTAL EQUIPMENT AND TEST CYLINDERS 15 

III . DATA ACQUISITION AND PROCESSING 19 

IV. GOVERNING PARAMETERS 20 

V. EXPERIMENTAL RESULTS 21 

A. CYLINDERS AT 90 DEGREE YAW ANGLE 21 

B, CYLINDERS AT 60 DEGREE YAW ANGLE 24 

C„ CYLINDERS AT 45 DEGREE YAW ANGLE 24 

Do COMPARISON OF THE DRAG AND INERTIA COEFFICIENTS 25 

E. EVALUATION OF MORISCN’S EQUATION FOR YAWED CYLINDERS 26 

Fo TRANSVERSE (LIFT) FORCES 27 

Go FOURIER ANALYSIS OF THE LIFT FORCE 27 

VI o CONCLUSIONS 28 

LIST OF REFERENCES 122 

INITIAL DISTRIBUTION LIST 123 



5 



LIST OF FIGURES 



1. A Schematic Drawing and Picture of the Water Tunnel 17 

2. versus K (1976 data) 30 

3. C m versus K (1976 data) 31 

4. versus K (1981 data) 32 

5. C m versus K (1981 data) 33 

6. versus K (1982 data) 34 

7. C m versus K (1982 data) 35 

8. versus K (1976, 1981, and 1982 data) 36 

9. C^ versus K (1976, 1981, and 1982 data) 37 

10. (a) First Mode of Vortex Shedding 38 

(b) Second Mode of Vortex Shedding 38 

11. versus K for a 6.5 inch Rough Cylinder (1976 data) 39 

12. C m versus K for a 6.5 inch rough Cylinder (1976 data) 40 

13. C^ versus K for a 6.5 inch rough Cylinder (1981 data) 41 

14. C m versus K for a 6.5 inch rough Cylinder (1981 data) 42 

15. C^ versus K for a 6.5 inch rough Cylinder (1982 data) 43 

16. C m versus K for a 6.5 inch rough Cylinder (1982 data) 44 

17. C^ versus K for a 6.5 inch rough Cylinder (1976, 1981, and 1982 

data) 45 

18. C m versus K for a 6.5 inch rough Cylinder (1976, 1981, and 1982 

data) 46 

19. C^ versus K for 8 = 3975, a = 90 Deg. , k/D = 0.00 47 

20. C m versus K for 3 = 3975, a = 90 Deg. , k/D =0.00 48 

21. C^ versus K for 3 = 4050, a = 90 Deg. , k/D = 0.01 49 

22. C m versus K for 8 = 4050, a = 90 Deg. , k/D = 0.01 50 

23. C^ versus K for 8 = 4057, a = 60 Deg. , k/D = 0.00 51 



6 



24 . C versus K for 0 = 4057 , a = 60 Deg . , k/D =0.00 52 

25. C^ versus K for 0 = 2271, a = 60 Deg. , k/D = 0.00 53 

26. C m versus K for 0 = 2271, a = 60 Deg. , k/D =0.00 54 

27. versus K for 0 = 1021, a = 60 Deg., k/D = 0.00 55 

28. versus K for 0 = 1021, a = 60 Deg. , k/D = 0.00 56 

29 . versus K for 0 = 4108 , a = 60 Deg . , k/D = 0.01 57 

30. C m versus K for 0 = 4108, a = 60 Deg. , k/D = 0.01 58 

31 . versus K for 0 = 2322 , a = 60 Deg . , k/D = 0.01 59 

32. C m versus K for 0 = 2322, a = 60 Deg. , k/D = 0.01 60 

33. C d versus K for 0 = 1042, a = 60 Deg., k/D = 0.01 61 

34 . C m versus K for 0 = 1042 , a = 60 Deg . , k/D =0.01 62 

35. C^ versus K for 0 = 4090, a = 45 Deg. , k/D = 0.00 63 

36. C m versus K for 0 = 4090, a = 45 Deg. , k/D = 0.00 64 

37. versus K for 0 = 2265, a = 45 Deg. , k/D = 0.00 65 

38. versus K for 0 = 2265, a = 45 Deg. , k/D = 0.00 66 

39. versus K for 0 = 1050, a = 45 Deg. , k/D = 0.00 67 

40 . C versus K for 0 = 1050 , a = 45 Deg . , k/D = 0.00 68 

41. versus K for 0 = 4150, a = 45 Deg. , k/D = 0.01 69 

42 . C m versus K for 0 = 4150 , a = 45 Deg . , k/D =0.01 70 

43. versus K for 0 = 2310, a = 45 Deg. , k/D = 0.01 71 

44. versus K for 0 = 2310, a = 45 Deg. , k/D =0.01 72 

45. C^ versus K for 0 = 1065, a = 45 Deg. , k/D = 0.01 73 

46 . C m versus K for 0 = 1065 , a = 45 Deg . , k/D =0.01 74 

47. Comparison of versus K for 0 = 4040, k/D = 0.00, a = 90, 

60, and 45 Deg 75 

48. Comparison of C^ versus K for 0 = 2345, k/D = 0.00, a = 90, 

60, and 45 Deg 76 



7 



49. 



Comparison of versus K for g = 1045, k/D = 0.00, 
a = 90, 60, and 45 Deg 77 

50. Comparison of C^ versus K for g = 4040, k/D =0.00, 

a = 90 , 60, and 45 Deg 78 

51. Comparison of C m versus K for g = 2345, k/D = 0.00, 

a = 90 , 60, and 45 Deg 79 

52. Comparison of C^ versus K for g = 1045, k/D =0.00, 

a = 90, 60, and 45 Deg 80 

53. Comparison of C^ versus K for g = 4100, k/D = 0.01, 

a = 90 , 60, and 45 Deg 81 

54. Comparison of C d versus K for g = 2325, k/D = 0.01, 

a = 90 , 60, and 45 Deg 82 

55. Comparison of C^ versus K for g = 1040, k/D = 0.01, 

a = 90 , 60, and 45 Deg 83 

56. Comparison of C^ versus K for g = 4100, k/D = 0»01, 

a = 90 , 60, and 45 Deg 84 

57. Comparison of C versus K for g = 2325, k/D = 0.01, 

a = 90, 60, and 45 Deg 85 

58. Comparison of C versus K for g = 1040, k/D = 0.01, 

a = 90 , 60, and 45 Deg 86 

59. Comparison of .Measured and Calculated Force (Two-Tern 
Morison Eauation) for K = 6.86, g = 4050, a = 90 Deg., 

k/D = 0.01 87 

60. Comparison of Measured and Calculated Force (Two-Term 
Morison Equation) for K = 10.18, g = 4050, a = 90 Deg. , 

k/D = 0.01 88 

61. Comparison of Measured and Calculated Force (Two-Term 
Morison Equation) for K = 11.88, g = 4050, a = 90 Deg., 

k/D = 0.01 89 

62. Comparison of Measured and Calculated Force (Two-Term 
Jforison Equation) for K = 16.15, g = 4050, a = 90 Deg., 

k/D = 0,01 90 



8 



63 . 



Comparison of Measured and Calculated Force (Two-Term 
Morison Equation) for K = 22,56, g = 4050, a = 90 Deg., 
k /D — 0,01 ,,,,«,000,.,00, 

64, Comparison of Measured and Calculated Force (Two-term 
Morison Equation) for K = 6,89, g = 4100, a = 60 Deg,, 
k/D =0,01 



9 * 9 * 9999*999 



• 009*0909 



65. Comparison of Measured and Calculated Force (Two-Term 
Morison Equation) for K = 16.83, g = 4100, a = 60 Deg,, 
k/D = 0,01 



90999 * 90 *<* 9999 * 9 * 9 * 999 > 



66, Comparison of Measured and Calculated Force (Two-Term 
Morison Equation) for K = 23.94, g = 4100, a = 60 Deg,, 
k/D = 0,01 



9999999*9 



99 9 * 999*9999 



67, Comparison of Measured and Calculated Force (Two-Term 
Morison Equation) for K = 7,90, g = 4150, a = 45 Deg,, 
k/D = 0,01 



9 9 9 9 



9 * 99*9 9 • 9 • • 



9 9 9 9 3 



68, Comparison of Measured and Calculated Force (Two-Term 
Morison Equation) for K = 10,05, g = 4150, a = 45 Deg,, 
k/D = 0,01 



999990 * 99999999 * 



69, Comparison of Measured and Calculated Force (Two-Term 
Morison Equation) for K = 12,50, g = 4150, a = 0.01 
k/D = 0,01 , 



90900099*0999 



9 9 9 * 



70, Comparison of Measured and Calculated Force (Two-Term 
Morison Equation) for K = 16,35, g = 4150, a = 45 Deg,, 
k/D = 0,01 



99990999999 * 99939 * 



71. Comparison of Measured and Calculated Force (Two-Term 
Morison Equation) for K = 24,04, g = 4150, a = 45 Deg., 
k/D = 0,01 , , 



999999*99 



9 * 999*9999 



72. Comparison of C- ^ ^ versus K for g = 4040, k/D = 0,00 
a = 90, 60, and 45 Deg, . . 



9 9 9 



999 * 9*99 



91 



92 



93 



94 



95 



96 



97 



98 



99 



100 



73, Comparison of C^^^ versus K for g = 4100, k/D = 0,01 

a = 90, 60, and 45 Deg, 101 



74, Normalized Lift Force Data for K = 8,38, 8 = 4090, 
a = 45 Deg. , k/D = 0,00 

75, Normalized Fourier Coefficients for K = 8.38, 8 = 4090, 
a = 45 Deg,, k/D =0,00 

76, Normalized Lift Force Data for K = 10,62, g = 4090, 
a = 45 Deg, , k/D =0,00 



. 102 
. 103 
, 104 



9 



77 . 



Normalized Fourier Coefficients for K = 10.62, 3 = 4090, 
a 45 Deg . , k / D 0 o 00 .oooo....oo.»ooo 



* 105 



78. Normalized Lift Force Data for K = 13.97, 3 = 4090, 

^ 45 Deg . , k /D 0 o 00 .oo..o.o«.oo.. 9 o. 

79. Normalized Fourier Coefficients for K = 13.79, 3 = 4090, 

a = 45 Deg. , k/D =0.00 ................. 

80. Normalized Lift Force Data for K = 17.63, 3 = 4090, 

ci 45 Deg . , k/D 0.00 ............o.... 

81. Normalized Fourier Coefficients for K = 17.63, 3 = 4090, 

ci = 45 Deg., k/D = 0.00 ................. 

82. Normalized Lift Force Data for K = 21.18, 3 = 4090, 

ci = 45 Deg. , k/D = 0.00 ................. 

83. Normalized Fourier Coefficients for K = 21.18, 3 = 4090, 

a = 45 Deg., k/D = 0.00 ................. 

84. Normalized Lift Force Data for K = 7.86, 3 = 4150, 

01 = 45 Deg., k/D = 0.01 ................. 

85. Normalized Fourier Coefficients for K = 7,86, 3 = 4150, 

ci = 45 Deg. , k/D = 0.01 ................. 

86 . Normalized Lift Force Data for K = 9.39 , 3 = 4150, 

a = 45 Deg., k/D = 0.01 ................. 

87. Normalized Fourier Coefficients for K = 9,39, 3 = 4150, 

ci 45 Deg. , k / L) — 0o01 . . . . . , . . . . , . . . . o « 

88 . Normalized Lift Force Data for K = 13.69, 3 = 4150, 

ci = 45 Deg., k/D = 0.01 ................. 

89. Normalized Fourier Coefficients for K = 13.69, 3 = 4150, 

01 = 45 Deg. , k/D = 0.01 

90. Normalized Lift Force Data for K = 17.57, 3 = 4150, 

01 = 45 Deg. , k/D = 0.01 ................. 

91. Normalized Fourier Coefficients for K = 17.57, 3 = 4150, 

ci = 45 Deg., k/D = 0.01 ................. 

92. Normalized Lift Force Data for K = 19,68, 3 = 4150, 

01 = 45 Deg, , k/D = 0.01 ................. 

93. Normalized Fourier Coefficients for K = 19.68, 3 = 4150, 

ci = 45 Deg., k/D = 0.01 ................. 



. 106 
, 107 
. 108 
. 109 
, 110 
. Ill 
, . 112 
, 113 
, . 114 
. . 115 
, , 116 
, . 117 
. . 118 
, . 119 
. . 120 
, , 121 



10 



TABLE OF SYMBOLS AND ABBREVIATIONS 
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I. INTRODUCTION 



The study of time dependent flow about yawed cylinders is a topic which 
presently has great interest, both theoretically and practically. Of the 
great variety of time-dependent flows about yawed cylinders, the wave motion 
and sinusoidally oscillating flow are of such importance as to demand inroe- 
diate investigation. For lack of adequate information concerning time- 
dependent flows, industry is often forced to adopt steady flow relations to 
time-dependent flow situations. This research proposes to investigate the 
forces acting upon yawed circular cylinders immersed in sinusoidally-oscil- 
lating planar flow and to examine the validity or limitations of the present 
methods of analysis. 

In 1950, Morison et al, [Refs, 1 and 2] introduced an equation for calcu- 
lating the in-line force acting on a vertical pile due to unbroken surface 
waves. For a cylinder of diameter D, the force per unit length is expressed 
as 

F = 0,5 pD C d U |U| + 0,25 TrpD 2 C m (DU/Dt) (1) 

where U represents the incident flow velocity; C^, the drag coefficient; and 

C , the inertia coefficient. The coefficients C and C, were considered time- 
rrr m d 

invariant and constant along the length of the cylinder, J.fcrison's paper was 
intended as a preliminary report with follow-on studies to be done on other 
structures in various wave actions. Morison et al. did not consider the 
contributions of transverse forces and vortex shedding in the calculation of 
the in-line forces. 

Attempts have been made to extend this rather simple relationship to 
far more coop lex situations such as combined wave and current flow, hydro- 
elastic oscillations, and yawed cylinders. Heideman et al, [Ref, 3] studied 
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the validity of Morison’ s equation using ocean data. Although Heideman et 
al, concluded that Morison ’ s equation was satisfactory for normal cylinders, 
the utilization of an equation unproven even under ideal flow situations 
[Ref, 2] to analyze ocean data so as to prove the validity of the equation 
is not a very meaningful exercise, (for additional details see Ref, 4), 

It is obvious that there is a great need for an idealization of the 
problem, or an experiment which is more manageable (e,g, , a sinusoidally 
oscillating planar flow). Only in this way can all the complex interactions 
be separately taken into account. 

Engineers faced with the problem of dealing with wave forces on yawed 
members and having no other recourse drew upon previous work with steady 
flow, Hoemer [Ref, 5] proposed the "independence principle", which stated 
that the normal pressure forces are independent of the tangential velocity 
for subcritical values of Re^, where Re n is the Reynolds number based upon 
the flow velocity normal to the cylinder. This principle allowed the decom- 
position of forces and velocities into normal and tangential components and 
the neglecting of the tangential conponents, Bursnall and Loft in [Ref, 6] 
found that the independence principle does not apply to the critical and 
transcritical flow regimes, Norton et al, [Ref, 7] found that the inde- 
pendence principle does apply to post -critical as well as subcritical flow, 
but not to the critical and transcritical regions in between. Thus, recent 
research has shown that the independence principle applies when the boundary 
layer is wholly laminar (Hoemer) or wholly turbulent (Norton), but its use 
in the critical and transcritical regions is uncertain. 

The designers of offshore structures were thus led to adopt the inde- 
pendence principle for the wave force calculations, and in doing so, to 
generalize the Morison equation. This has been accepted practice in industry 
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for lack of a better relation. It can be asked what would be anticipated 

for oscillating flow or waves on the basis of what is known for steady flow. 

The instantaneous Reynolds number in such a flow will vary from -Re_ to 

max 

+ R e nax during a complete flow cycle. It could be postulated that the bound- 
ary layer would, at times, be fully laminar; at other times fully turbulent; 
and the rest of the time be in transition. In light of this, it is rather 
doubtful that the independence principle applies at all to oscillating or 
wave flow. 

The study of forces acting on yawed cylinders can be carried out 
either by oscillating a yawed cylinder in a tank, or using small amplitude 
waves in a laboratory channel, oscillating the flow about a fixed cylinder, 
or by using ocean data if available. Oscillating the cylinder has proven 
to be impractical [Ref, 8] because of difficulties of accounting for the 
inertial force acting on the body, producing repeatable oscillations, 
vibrations in the system, and difficulties in measuring the in-line and 
transverse forces simultaneously. Waves have relatively more complex 
flow kinematics due to the orbital motion of particles and the decay of 
wave amplitude with depth. Ocean data are not available for such a study. 
Oscillation of the flow past the cylinder has proved to be the best method 
[Refs, 9 and 10], 

The primary objective of this investigation was to study the forces 
exerted by a sinusoidally oscillating planar flow on yawed circular cylinders 
to determine whether the independence principle is applicable or not. If so, 
the force transfer coefficients calculated by Fourier analysis using the 
normal velocity component should reduce identically to the normal cylinder 
case at corresponding values of K, Re, and k/D, If the independence prin- 
cipal does not apply, it is desired to determine what the coefficients are 



14 



as functions of yaw angle, roughness ratio k/D, Re, and K. It would also 
be necessary to determine how well Morison's equation works with the new 
force coefficients. 

With the foregoing objectives in mind a detailed investigation has 
been undertaken using smooth and rough cylinders of nominal 6 inch, 4.5 
inch, and 3 inch diameters at yaw angles of 45 degrees, 60 degrees, and 90 
degrees (as measured between the ambient flow direction and the cylinder 
axis). The relative sand-roughness for all rough cylinders was k/D = 1/100, 



II. EXPERIMENTAL EQUIPMENT AND TEST CYLINDERS 

The equipment consists of a large U-shaped oscillating flow tunnel. 

It was first constructed in 1975. Since then the length of the tunnel has 
been increased from 30 feet to 35 feet and its height from 16 feet to 22 
feet. The cross-section of the 35 ft long horizontal test section has been 
increased from 3 ft by 3 ft to 3 ft by 4,7 ft. Furthermore, the oscillation 
mechanism has been conpletely modified so that mono-harmonic oscillations 
can be generated and maintained indefinitely at the desired amplitude. For 
this purpose the output of a 2 Hp fan was connected to the top of one of 
the legs of the tunnel with a large pipe (D = 3 ft), A small butterfly 
valve, placed in a special housing, outside the tunnel, between the top of 
the tunnel and the supply line, oscillated continuously and sinusoidally at 
a frequency equal to the natural frequency of the oscillations of water in 
the tunnel. The oscillation of the valve was perfectly synchronized with 
that of the flow through the use of an electronic feedback control system, 
coupled to a DC motor oscillating the valve plate. The circuit maintained 
the period of oscillations of the valve within 0,0005 seconds. The amplitude 
of the oscillations was varied by constricting or enlarging an orifice at the 
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exit of the fan* The flow oscillated at a given anplitude as long as 
desired, (see Fig. 1 for a picture and schematic drawing of the tunnel ) . 

The velocity in the tunnel has been determined through the use of 
a capacitance wire, hotfilm anemometer, perforated ball, magnetic flow 
meter, an accelerometer (which measured the instantaneous acceleration of 
flow in the test section), and by visual measurement of the water level 
at its heighest and lowest points in the legs of the tunnel. It is safe 
to state that the velocity could not have been measured mere accurately. 
The only other means by which the velocity could have been measured was 
the use of a laser device. In view of its cost and in view of the fact 
that the other means of measurement yielded the ambient velocity within 
2 percent of each other it was decided to forsake the laser system. 

Mounting pads for the strain gage housings were placed on opposite 
sides of the tunnel in such a manner so as to accommodate the yawed cyl- 
inders. A separate mounting had to be installed for each angle of yaw. 
The mounting pad on one side of the tunnel was fitted with an adjustable 
slide so that the pad could be moved horizontally small distances to 
exactly match the cylinder length. The force transducers and housings 
were unchanged from descriptions in [Refs. 9 and 10]. 

The cylinder ends, cut exactly at the desired angle of yaw, were 
parallel to the tunnel walls with a 1/32 inch gap at each end. This gap 
was filled with a soft foamy material glued to the ends of the cylinder. 
Ball bearings were mounted in the ends of the cylinder with the outer 
bearing faces flush with the face of the cylinder. Careful calibration 
and extensive testing proved that this mounting system allowed accurate 
and repeatable recording of the normal and transverse forces upon the 
test cylinder. 
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Fig. 1 A schematic drawing and picture of the water tunnel 
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The calibration of the cylinder was first conducted in a vertical 
direction, since yaw would not affect the total force in this direction. 
Weights were hung from the center of the cylinder, and calibration factors 
were obtained by converting the sum of the signals of two force trans- 
ducers to either pounds/nm or pounds/volt. The strain gages were found 
to be exactly linear throughout the range of expected forces . Ideally, 
the horizontal forces normal to the cylinder should be related to the 
streamwise force by the sine of the yaw angle a (as measured between the 
cylinder axis and the ambient flow direction) . A system of supports and 
pulleys allowed applying horizontal force at the center of the cylinder and 
normal to it. The sum of the outputs of the two transducers was recorded 
on the strip chart recorder, establishing the necessary relationship between 
the normal in-line force F and the total electrical output of the gages. 

The same gages were used to measure the lift force or the in-line force 
by rotating the gages 90 degrees. It is easy to show that the calibration 
factor for the lift force (here in the vertical direction, up or down), 
expressed in terms of pounds or millimeter deflection, is sina times the 
calibration factor for the in-line force. This proved to be true experi- 
mentally also and demonstrated independently the validity of the measure- 
ment technique. 

The following cylinder sizes and angles were tested: 

90 degrees, 6.5 in. and 6 in., anooth and rough; 

60 degrees, 6.0 in., 4.45 in., and 3 in., smooth and rough; 

45 degrees, 6 in., 4.45 in., 3 in., smooth and rough. 

Relative sand roughness of k/D = 1/100 was used for all rough cylinders. 
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III. DATA ACQUISITION AND PROCESSING 



The data were acquired using an HP-3052A automatic data acquisition 
system. The system consisted of an HP- 349 5 A 20 channel scanner and an 
HP-3437A system voltmeter, both controlled by an HP-9845B desk top com- 
puter. The in-line and transverse force signals, originating as a 
voltage fran the strain gages, was then amplified and sent to the strip 
chart recorder. This amplified signal was also sent to the scanner as 
one channel of data. The flow amplitude signal, originating from a 
differential pressure transducer, was similarly amplified and sent to 
both the strip chart recorder and the scanner as another channel of data. 
The strip chart data was maintained as a visual record of the flow ampli- 
tude and force. The scanner switched between the two channels at desig- 
nated intervals when triggered by the voltmeter. The voltmeter read each 
one of these voltage values and transferred them to the computer for 
immediate calculation of the governing parameters or for further trans- 
fer to floppy disc storage to allow for analysis at a later time. 

At least six cycles of in-line force data, digitized at 0.5 degree 
intervals, were acquired with the HP-3052A - HP-9845B system. The data 
were then averaged to give one cycle from which the governing parameters 
were calculated. Ten cycles of transverse force data, digitized at one 
degree intervals, were acquired with the same system. These data were not 
averaged before calculating the governing parameters. 

The use of the data acquisition system proved to be highly beneficial 
for a number of reasons. The calculated parameters were available on a 
real time basis, allowing for detailed study of anomalies in the force 
patterns. Taking data over a larger number of cycles removed more of the 
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small randomness of the flow forces. This was particularly valuable for 
analyzing the transverse force data. 



IV. GOVERNING PARAMETERS 



Data reduction for the in-line forces is based on Morison's equation 
using the normal component of velocity. The derivation of the governing 
parameters is described in detail in [Refs. 2, 9-11]. Here only a brief 
description of the parameters will be given. 

The force per unit length of a cylinder of diameter D is represented 
by Morison's equation as 



C = — = — C sin0 - C , I cos9 | cos0 
T Z K. m d 1 



( 2 ) 



pDU 



m 



For an oscillating flow represented by U = -U^cos© with 0 = 27rt/T, 
where U m represents the maximum flow velocity, the Fourier averages of C m 



and are given by 



2tt 



F sin0 



C m = 



= 2K f r rn 

it 3 J p D l 



U 2 



d0 



(3) 



ra 



and 
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q r 

C = _ 1. I JL 
d 4 / n 

J p D 



F„cos0 



d0 



(4) 



pDD, 



0 



m 



where F represents the measured force, 
m 

The transverse force data are represented by the normalized rms 
value of the measured lift force (C ^ „^ ) . The lift data were also analyzed 
using a Fourier transform to obtain frequency information. The first 15 
Fourier coefficients, with corresponding magnitude and phase angles, were 
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calculated. Finally, C^, C , and C 1 rmc; were plotted in terms of the 

2 

Keulegan-Carpenter number K, the 'frequency parameter' 3 = Re/K = D /vT, 
the relative roughness k/D, and the yaw angle a. 



V. EXPERIMENTAL RESULTS 

A. CYLINDERS AT 90 DEGREE YAW ANGLE 

Attempts to achieve as high Reynolds numbers as possible in conducting 
wind-tunnel and water-tunnel experiments invariably give rise to wall- 
interference effects which, of course, influence whatever measursnents are 
made. There are several blockage correction formulas for steady flows 
which might be used so that the wall-interference effects on the calculated 
force coefficients might be minimized. Unfortunately, none of these formulas 
could be used in the present study for no one has danonst rated that the 
blockage effects in oscillatory flows are identical to those experienced in 
steady flows [Ref. 9]. 

In view of the foregoing it was decided to repeat the experiments 
reported in 1976 [Refs. 9 and 10] with a 6.5 inch smooth and sand-roughened 
cylinder (k/D = 1/100). Originally the experiments were conducted in the 
first version of the U-shapped water tunnel which had a length of 30 ft, a 
height of 16 ft, and a cross-section of 3 ft by 3 ft. Subsequently, the 
length of the tunnel has been increased frcm 30 ft to 35 ft and its height 
fran 16 ft to 22 ft. The cross-section of the 35 ft long horizontal test 
section has been increased from 3 ft by 3 ft to 3 ft by 4.7 ft, as noted 
earlier. Finally, the data acquisition system and other refinements have 
been incorporated into the system. Experiments with the 6.5-inch cylinder 
have been repeated following each change in the tunnel. 
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The drag and inertia coefficients obtained in 1976 with a D = 6.5 in. 
smooth cylinder are shown in Figs. 2 and 3. The data obtained with the 
same cylinder in the larger tunnel in 1981 and in 1982 are shown in Figs, 

4 through 7. Finally, the entire data for and C m are shown in Figs. 8 
and 9. Evidently, the data obtained with the enlarged tunnel agree 
extremely well with those obtained in 1976 [Ref. 9]. No clearer proof 
can be presented than these figures in eliminating the possibility of any 
blockage effect on the data reported previously by Sarpkaya [Refs. 9 and 
10], particularly at large values of K and Re. 

Figures 8 and 9 show that the only difference between the 1976 data and 
those obtained in 1981 and 1982 is the occurrence of two distinct modes of 
vortex shedding for a given K in the region of K values from 9 to 16, The 
observations of the in-line and transverse-force traces over several hundred 
cycles of flow oscillation in the said region of K values (the drag-inertia 
dominated regime) have revealed that the dominant mode of flow is as shown 
in Fig. 10a. For this mode, the in-line force trace does not exhibit large 
second order oscillations and the transverse force is both small and irregular. 
In other words, at such large Reynolds numbers (Re varied from about 47,000 
to 85,000 as K varied from 9 to 16) the spanwise coherence of vortices is 
not perfect. The in-line force corresponding to this dominant mode yields 
drag and inertia coefficients identical to those shown in Figs. 2 and 3 
(the 1976 data). From time to time, however, the spanwise coherence of the 
vortices improves dramatically and the the in-line and transverse forces 
change to those shown in Fig. 10b, In this second mode of flow, the in-line 
force exhibits large oscillations and the magnitude of the transverse force 
increases dramatically. 'Ibis second mode gives rise to the larger drag and 
amller inertia coefficients shown in Figs. 8 and 9 in the region 9 < K < 16. 
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Thus, it is clear that one can obtain two different drag coefficients and 
hence two different inertia coefficients for the same K (in the drag- 
inertia dominated regime), depending on the flow mode. The occurrence of 
the second mode in the drag-inertia dominated regime for Reynolds numbers 
smaller than about 30,000 is quite cornnon [Ref. 9], As the Reynolds number 
increases, however, the first mode becomes the dominant mode. Then the 
variation of the transverse force with time becomes a non-st at ionary random 
process. 

It is clear from the foregoing that unique values of C d , C m , and 
do not exist for certain values of K and Re. It is also clear that the 
orbital motion of the fluid particles and the omnidirectionality of the 
waves and currents in the ocean environment give rise to the first flow 
mode and result in smaller drag and lift coefficients. In laboratory 
experiments the coherence length emerges as an important parameter, 
particularly in the drag- inertia dominated regime. It is thus possible to 
choose and L/D radio for which the flow will be in the first or second mode 
for a given combination of K and Re in the drag- inertia dominated regime. 

The drag and inertia coefficients obtained with a D = 6.5 inch rough 
cylinder (k/D = 1/100) in 1976 are shown in Figs. 11 and 12. The data 
obtained with the same cylinder in 1981 and 1982 are shown in Figs, 13 
through 16. Finally, a composite plot of the entire data is shown in Figs, 

17 and 18. Evidently, all three sets of the data agree extremely well 
and show that tunnel blockage did not play any role on the data reported 
previously by Sarpkaya [Refs. 9 and 10] for smooth as well as rough cylinders. 

The observations of the in-line and transverse force traces for the 
rough cylinder have shown that only the second flow mode occurred in the 
drag- inertia dominated regime. Evidently, roughness increases the coherence 
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length, giving rise to larger drag coefficients and smaller inerti coeffi- 
cients in the drag-inertia dominated regime. Additional work is needed 
to delineate more clearly the role played by the coherence length and 
surface roughness in the variation of the force-transfer coefficients. 

The drag and inertia coefficients obtained with a D = 6 inch nominal 
diameter smooth and rough cylinder are shown in Figs, 19 through 22, The 
results are quite similar to those presented earlier for the D = 6,5 inch 
smooth and rough cylinder, 

B, CYLINDERS AT 60 DECREE YAW ANGLE 

The drag and inertia coefficients obtained with 6 inch, 4.5 inch, and 
3 inch (nominal diameter) smooth and rough cylinders (k/D = 1/100) are shown 
in Figs, 23 through 34 as a function of K (K = U T/D). Two facts are 
inmediately apparent. First, the drag coefficient does not exhibit multiple 
values in the drag- inertia dominated regime. In general, the drag coefficient 
increases to a peak and then decreases gradually with increasing K and Re, 

In fact, the observations of flow about the cylinder with tracer particles 
have shown that the vortices are far from coherent. Secondly, the inertia 
coefficient is larger than 2,0 for small values of K, giving the first 
indication that the independence principle is not probably valid, 

C, CYLINDERS AT 45 DEGREE YAW ANGLE 

The drag and inertia coefficients obtained with 6 inch, 4.5 inch, and 
3 inch (nominal diameter) smooth and rough cylinders (k/D = 1/100) are shown 
in Figs. 35 through 46 as a function of K. Aside from a relatively larger 
scatter, the overall characteristics of the data are quite similar to those 
for the 60 degree yaw angle. 
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D. COMPARISON OF THE DRAG AND INERTIA COEFFICIENTS 



The drag coefficient for three smooth cylinders at yaw angles of 90, 

60, and 45 degree are shown in Figs. 47 through 49. Evidently, had the 
independence principle been valid, the three sets of data corresponding to 
three yaw angles would have collapsed into a single set, thus rendering C^ 
versus K relationship independent of the yaw angle (at least within the 
range of the test parameters and experimental errors) . 

The inertia coefficient for the three smooth cylinders at yaw angles 
of 90, 60, and 45 degree are shown in Figs. 50 through 52. Evidently, C m 
varies significantly with the yaw angle and the independence principle is 
not valid within the range of test parameters. The inertia coefficients 
based on the normal component of acceleration are significantly larger than 
those for the normal cylinder (yaw angle = 90 degrees) at the corresponding 
Reynolds numbers and Keulegan-Carpenter numbers. Finally, the inertia 
coefficient for the yawed cylinders does not exhibit the usual "inertia 
crisis" associated with the normal cylinder. In fact, the inertia coeffi- 
cient for the yawed cylinders varies gradually with K and remains 25 percent 
to 50 percent larger than that for the normal cylinder. Thus, it is clear 
frcm the foregoing that the normal force acting on a yawed cylinder is 
significantly underestimated through the use of Morison's equation, inde- 
pendence principle, and the drag and inertia coefficients appropriate to the 
normal cylinder. 

The drag and inertia coefficients for the three rough cylinders at yaw 
angles of 90, 60, and 45 degree are shown in Figs. 53 through 58. Evidently, 
roughness plays an important and, oddly enough, a unifying role on the drag 
coefficient. This is somewhat anticipated on the basis of the fact that 
roughness precipitates earlier transition and greater incoherence in the 
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flow about a yawed cylinder, thereby creating conditions more favorable 
to the independence principle. Figures 53 through 55 show that the inde- 
pendence principle for rough cylinders is almost valid with the exception 
of the drag- inertia dominated regime. However, the inertia coefficient 
for rough cylinders depends strongly on the yaw angle (see Figs. 56 through 
58). Thus, it appears that the independence principle or the cosine law, 
as it is sometimes called, is a gross simplification of the behavior of 
flow in the near wake. One may, therefore, conclude that the Fourier- 
averaged drag and inertia coefficients, based on Morison's equation, are 
unique for each angle of yaw, Reynolds number, Keulegan -Carpenter number, 
and the relative roughness. 

E. EVALUATION OF MORISON'S EQUATION FOR YAWED CYLINDERS 

Since the independence principle cannot be applied to yawed cylinders 
in sinusoidally oscillating planar flow, it is necessary to determine if the 
Fourier-averaged coefficients calculated frcm the normal forces will predict 
the forces exerted on the cylinder when Morison's equation is used together 
with the experimentally determined C^ and C^ values, appropriate to each yaw 
angle. In general it was found that the Morison equation predicts the forces 
on the yawed cylinder with the same degree of accuracy as the 90 degree case. 
Figures 59 through 63 show comparisons of the measured and calculated forces 
for the 90 degree 6-inch rough cylinder at K values of 6.86, 10.18, 11.88, 
16.15, and 22.56. Examples for the 60 degree case are shown in Figs. 64 
through 66 for K values of 6.89, 16.83, and 23.94. Finally, Figs. 67 through 
71 show the comparison of the measured and calculated forces for the 45-degree 
case for K values of 7.90, 10.05, 12.50, 16.35, and 24.04. It is clear from 
these comparisons that when C and are known for a given angle of yaw, 
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Morison ' s equation will predict the normal forces on cylinders at angles 
of yaw as well as the normal cylinder case. 

F. TRANSVERSE (LIFT) FORCES 

The lift force data are presented in terms of the normalized rms value 
of the measured lift force. Figures 72 and 73 show the lift coefficients 
for the 6- inch smooth and rough cylinder, respectively, for the three yaw 
angles. The lift coefficient for the 60 degree and 45 degree yaw angles 
is considerably smaller than that for the normal smooth cylinder. Further- 
more, the data for the 60-degree and the 45-degree cases exhibit great deal 
of scatter. This is primarily due to the fact that the coherence length is 
considerably smaller and the flow exhibits numerous modes, leading to a 
non-st at ionary random process as far as the vortex shedding is concerned. 

For the rough cylinder, the lift coefficients are relatively larger 
for all values of K. However, the lift coefficient for the 60-degree and 
45-degree cases is about half of that for the normal cylinder for values 
of K larger than about 15. The data for all other smooth and rough cylinders 
showed similar trends and will not be discussed here further. 

G. FOURIER ANALYSIS OF THE LIFT FORCE 

Calculation of the first fifteen Fourier coefficients, magnitudes, and 
phase angles of the lift force only reinforced the fact that the lift force 
is a highly random process. Except for low values of K (where second harmonic 
was dominant) there was no clear dominance of any particular harmonic. In 
fact , it showed that a number of harmonics were present . Neither was any 
trend visible in the phase angle variation. Figures 74 through 93 show sample 
data plots, and harmonics for the 45-degree 6-inch rough and smooth cylinders 
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at various K values. The existence of a large number of harmonics in the 

lift force accounts for most of the scatter in CL 

lrms 

VI. CONCLUSIONS 

The forces acting on yawed smooth and rough circular cylinders in 
sinusoidally oscillating planar flow have been investigated extensively 
and the following conclusions have been reached: 

1, The independence principle does not apply over the range of K 
and Re values investigated. The drag and inertia coefficients for the 
45-degree and 60-degree smooth and rough yawed cylinders differ signifi- 
cantly from those for the 90-degree normal cylinder, 

2, The Fourier-averaged drag and inertia coefficients, based on 
Morison’s equation, and the rms value of the lift coefficient are unique 
for each yaw angle, Reynolds number, Keulegan-Carpenter number, and the 
relative roughness, 

3, When the drag and inertia coefficients are calculated from the 
measured force for the yawed cylinder, in a manner similar to that used 
for the 90-degree cylinder, Mori son ' s equation predicts the measured force 
with the same accuracy as that for the 90-degree cylinder, 

4, The normal force acting on a smooth or rough yawed cylinder is 
significantly underestimated through the use of the independence principle 
and the drag and inertia coefficients appropriate to the normal cylinder, 

5, Extensive flow visualization about normal and yawed cylinders has 
shown that the cylinder inclination significantly decreases the spanwise 
coherence and the kinematics of flow about yawed cylinders are considerably 
more complex than those for a normal cylinder. 
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6, The lift force is a non-st at ionary random process and contains a 
large number of harmonics „ In general, for yawed cylinders is 

considerably smaller than that for the normal cylinder, 

7„ The data presented herein should form the basis of future 
calculations for the forces acting on yawed cylinders. Additional data 
at higher Reynolds numbers and the understanding of the role played by 
the coherence length will add significantly to the quantitative and 
qualitative -understanding of flow about yawed cylinders. 
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Figure 2. C, versus K (1976 data) 
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Figure 3. C versus K (1976 data) 
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Figure 4. C , versus K (1981 data) 
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Figure 5. C versus K (1981 data) 




<V1 CVJ — * — * * - 



34 



Figure 6. C, versus K (1982 data) 
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Figure 7. C versus K (1982 data) 
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Figure 8. C, versus K (1976, 1981, and 1982 data) 
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Figure 9. C versus K (1976, 1981, and 1982 data) 
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Figure 11. C , versus K for a 6.5 inch rough cylinder (1976 data) 
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Figure 12. C versus K for a 6.5 inch rough cylinder (1976 data) 




«l N - - — 



41 



Figure 13. C, versus K for a 6.5 inch rough cylinder (1981 data) 
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Figure 14. C versus K for a 6.5 inch rough cylinder (1981 data) 
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Figure 15. C, versus K for a 6.5 inch rough cylinder (1982 data) 
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Figure 16. C versus K for a 6.5 inch rough cylinder (1982 data) 
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Figure 17. C versus K for a 6.5 inch rough cylinder (1976, 1981, and 1982 data) 
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Figure 18. C versus K for a 6.5 inch rough cylinder (1976, 1981, and 1982 data) 
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Figure 19. C. Versus K for /9= 3975, (X = 90 Deg., k/D = 0.00 
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Figure 20. C Versus K for p= 3975, 90 Deg., k/D = 0.00 
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Figure 21. C, Versus K for 8 = 4050, <X= 90 Deg., k/D 
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Figure 22. C Versus K for p= 4050, 0C= 90 Dec., k/D 
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Figure 23. C, Versus K for 3 = 4057, a = 60 Deg., k/D = 0.00. 
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Figure 25* C, Versus K for 3 — 2271, a - 60 Deg., k/D — 0.00. 
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Figure 27 . C, Versus K for g = 1021, a = 60 Deg., k/D 
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Figure 31. C , Versus K for 3 = 2322, a = 60 Deg., k/D = 0.01. 
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Figure 32. C Versus K for 3 = 2322, a = 60 Deg., k/D = 0.01. 
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Figure 33. C, Versus K for 3 = 1042, a = 60 Deg., k/D = 0.01. 
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Figure 37. C, Versus K for /}= 2265, 0C= 45 Deg., k/D = 0.00 
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Figure 39. C, Versus K for 18 = 1050, 0t= 45 Deg., k/D = 0.00 
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Figure 40. C Versus K for /3= 1050, 01 = 45 Deg., k/D = 0.00 
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Figure 41. C, Versus K for p= 4150, 0t= 45 Deg., k/D = 0.01 
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Figure 42. C Versus K for p= 4150, U= 45 Deg., k/D = 0.01 
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Figure 43. C, Versus K for f&= 2310, 0t= 45 Deg . , k/D = 0.01 
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Figure 44. C Versus K for /8= 2310, & = 45 Deg., k/D =0.01 
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Figure 45. C, Versus K for D= 1065, 0t= 45 Deg., k/D = 0.01 
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Figure 46. C Versus K for £= 1065, 0L= 45 Deg., k/D 
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Figure 49. Comparison of C-, Versus K for 3 = 1045, k/D = 0.00, 
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Figure 51. Comparison of C Versus K for 3 = 2345, k/D = 0.00, 
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Figure 52. Comparison of Versus K for 3 = 1045, k/D = 0.00 
a = 90, 60, and 45 Deg. 
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Figure 53. Comparison of C, Versus K for 3 = 4100, k/D = 0.01, 
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Figure 55. Comparison of Versus K for g = 1040, k/D = 0.01 
a = 90, 60, and 45 Deg. 
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Figure 57. Comparison of Versus K for £ = 2325, k/D = 0.01, 
a = 90, 60, and 45 Deg. 
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Figure 59. Comparison of Measured and Calculated Force 
(Two-Term Morison Equation) for K = 6.86, 
jQ= 4050, CX= 90 Deg., k/D = 0.01 
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Figure 60. Comparison of Measured and Calculated Force 
(Two-Term Morison Equation) for K = 10.18, 
/ 3 = 4050, 0t= 90 Deg., k/D = 0.01 
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Figure 61. Comparison of Measured and Calculated Force 
(Two-Term Morison Equation) for K = 11.88, 
/3= 4050, OC= 90 Deg., k/D = 0.01 
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Figure 62. Comparison of Measured and Calculated Force 
(Two-Term Morison Equation) for K = 16.15, 
/3= 4050, OL= 90 Deg., k/D = 0.01 
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Figure 63 . Comparison of Measured and Calculated Force 
(Two-Term Morison Equation) for K = 22.56, 
4050, 01 = 90 Deg., k/D = 0.01 
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Figure 64 . Comparison of Measured and Calculated Force 
(Two-Term Morison Equation) for K = 6.89, 

/Q = 4100, Oi= 60 Deg., k/D - 0.01 
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Figure 65. Comparison of Measured and Calculated Force 
(Two-Term Mori son Equation) for K = 16.83, 
/3 = 4100, OC= 60 Deq., k/D = 0.01 
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Figure 66 . Comparison of Measured and Calculated Force 
(Two-Term Morison Equation) for K - 23.94, 
/ 3 = 4100, OL= 60 Deg., k/D = 0.01 
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Figure 69. Comparison of Measured and Calculated Force 
(Two-Term Morison Equation) for K = 12.50, 

B = 4350, OL= 45 Deg., k/D = 0.01 



FORCE Vs RNGLE (deg) 



09S 




« « 
(73 (XI 



Q 



• • 

oj cn 

i i 



98 



Figure 70. Comparison of Measured and Calculated Force 
(Two-Term Morison Equation) for K = 16.35, 

/3 = 4150, 0i= 45 Deg., k/D = 0.01 
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Figure 71 . Comparison of Measured and Calculated Force 
(Two-Term Morison Equation) for K = 24.04, 
4150, Oi= 45 Deg., k/D = 0.01 
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Figure 72. Comparison of Versus K for 0 = 4040, k/D = 0.00, 
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Figure 73. Comparison of Versus K for 3 = 4100, k/D = 0.01 
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Figure 74. Normalized Lift Force Data for K = 8.38 
/9= 4090, OC= 45 Deg., k/D = 0.00 
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Beta = 4075 
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Figure 81. Normalized Fourier Coefficients for K = 17.63 
/9= 4090, 45 Deq., k/D = 0.00 
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Figure 83. Normalized Fourier Coefficients for K = 21.18, 
>3 = 4090, oL= 45 Deg., k/D = 0.00 
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Figure 84. Normalized Lift Force Data for 

4150, &L = 45 Deg., k/D = 0 
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Figure 86. Normalized Lift Force Data for K = 9:39 
&= 4150, 04 = 45 Deg., k/D = 0.01 
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Figure 87. Normalized Fourier Coefficients for 

= 4150, Oi= 45 Deg., k/D = 0.01 
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Figure 88. Normalized Lift Force Data for K = 12.69 
/$= 4150, 6 L= 45 Deg., k/D = 0.01 
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Figure 90. Normalized Lift Force Data for K = 17.57 
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Beta = 4162 
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Figure 93. Normalized Fourier Coefficients for K = 19.68, 
/S= 4150, 45 Deg., k/D = 0.01 
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